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Abstract 



I Using the Nahm transform we investigate doubly periodic charge one SU{2) in- 

• stantons with radial symmetry. Two special points where the Nahm zero modes have 

, softer singularities are identified as the locations of instanton core constituents. For a 

I square torus this constituent picture is closely reflected in the action density. In rect- 

angular tori with large aspect ratios the cores merge to form monopole-like objects. 
Oh. For particular values of the parameters the torus can be cut in half yielding two copies 

(— I I of a twisted charge ^ instanton. These findings are illustrated with plots of the action 

density within a two-dimensional slice containing the constituents. 

X. 

■ 1 Introduction 

Monopoles are at the heart of various kinds of electric- magnetic duality in extended super 
Yang-Mills theory, which have led to an impressive body of quantitative results. The role 
of monopoles and other topological constructs in non-supersymmetric Yang-Mills theory is 
less clear. Using such objects, particularly monopoles and vortices, models of confinement 
and chiral symmetry breaking have been advanced. A weakness of these models is that clas- 
sical Yang-Mills does not possess monopole solutions, rather they have instanton solutions. 
However, in recent years it has become clear that periodic instantons have a substructure 
that can be identified with monopoles and other extended objects. In particular, charge 
one SU (N) calorons (finite temperature instantons) can be viewed as bound states of 
monopole constituents P3I2]- Provided the 'size' of the instanton is larger than the period 
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(or inverse temperature) tubes of action density form. These can be identified with the 
worldhnes of the constituent monopoles. It would be interesting to extend these results to 
higher charge sectors and higher tori, i.e. more than one period so that we are considering 
instantons on T" x n > 1. Indeed, lattice based numerical studies indicate that dou- 

bly periodic instantons (n = 2) have remarkable properties including fractionally charged 
solutions, vortex-like objects and an exponentially decaying action density None of these 
attributes are shared by the calorons. Nevertheless, it would be desirable to have a unified 
description of calorons and multiply-periodic instantons and their constituents. 

The calorons are best understood within a formalism due to Nahm |31 13] . This is an 
extension of the ADHM construction. Nahm's approach was generalised to other four man- 
ifolds, notably the four-torus, 0. Here it becomes a duality, mapping U{N) charge k 

~ 4 

instantons on the four torus to U{k) charge N instantons on the dual torus, T (defined by 
inverting the four periods). The caloron construction and even ADHM can be understood as 
a limiting case of the four torus duality. One therefore expects that doubly periodic instan- 
tons can be treated in a similar fashion. Jardim [Zj has used the Nahm transform to discuss 
the charge one case. In [H] doubly periodic charge one instantons with radial symmetry 
in the non-compact directions were considered. Under the Nahm transformation they 
are mapped to abelian potentials on the dual torus T . Starting with these rather simple 
Nahm potentials the original instantons can be 'recovered' via the inverse Nahm transform. 
Technically, this involves solving certain Weyl-Dirac equations; for each x G x one has 

~ 2 

a different Weyl equation on T . The Weyl zero modes were determined explicitly for a two 
dimensional subspace (corresponding to the origin of M^). 

In a recent letter PJ we used these zero modes to compute the action density, analytically 
and numerically, within the subspace. We found that there are points in x where the 
solution of the Weyl equations have softer singularities. These are obvious candidates for 
constituent locations. We interpreted these constituents as overlapping instanton cores. For 
square tori our constituent picture correlates closely to the action density. Some of the charge 
one solutions we considered can be cut in half yielding two copies of a charge | instanton. 
This procedure is only consistent if the gauge potential has specific asymptotic boundary 
conditions which are analogous to a four torus twist. 

In this paper we obtain further information about doubly periodic instantons. In partic- 
ular, action density plots for rectangular tori are presented. These clearly show how doubly 
periodic instantons interpolate between fractionally charged lumps characteristic of four 
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torus solutions and monopole like tubes of actions density typical of calorons. Furthermore, 
we present an analytical discussion of the exponential decay of the action density. The paper 
is organised as follows. In section 2 the basics of the Nahm transform are briefly reviewed, 
and its application to doubly periodic instantons is discussed. Sections 3 to 5 deal with the 
Nahm potentials, zero modes and constituents corresponding to charge one instantons with 
gauge group SU{2). Section 6 describes some technical results regarding Green's functions 
that crop up in field strength calculations. Our numerical action density computations are 
detailed in section 7. Section 8 considers charge ^ solutions and the asymptotics of the gauge 
field and field strength are discussed in section 9. 

2 Nahm Transform for Doubly Periodic Instantons 

Before we specialise to x let us briefly recall how the Nahm transform is formulated 
on the four torus jO] (for a string-theoretic approach see [E]). Start with a self-dual anti- 
hermitian SU{N) potential, Afj_{x), on a euclidean four-torus with topological charge k. 
A gauge field on is understood to be an potential which is periodic (modulo gauge 
transformations) with respect to — x^+L^ (/i = 0, 1, 2, 3), the being the four periods 
of the torus. Self-duahty is the requirement that F^^, = ^e^^a/s^a/s, and we have taken the 
convention 60123 = —1- The next step is to turn the SU{N) instanton into a U{N) instanton 

by adding a constant U{1) potential, A^{x) — > ^4^(2;) — iz^, where the are real numbers. 

~ 4 

We can regard the 2;^ as coordinates of the dual torus , T , since the shifts Zf^ ^ z^ + 2n / 
can be effected via periodic U{1) gauge transformations. 
Now consider the U (N) Weyl operator 

D^iA) = cr^D'^iA), Di^iA) = + A^(x) - iz^ (1) 

with cr^ = (1, iri, ?r2, irs) where the Xj are Pauli matrices. Provided certain mathematical 
technicalities are met Dl{A) = —ajj^D'^{A) has k square-integrable zero modes iIj^{x] z) with 
i = 1,2, k. The Nahm potential is defined as 

A^^iz) = I^J^x4j'\x;z)-^^^{x;z). (2) 

Here the zero modes are taken to be orthonormal. Remarkably, A{z) is a U{k) instanton on 
the dual torus with topological charge A^. The associated field strength can be written as 
follows 

^;i(^)=/" d'^ f d^x'^'\x-z)a^{DlD,)-\x,x')al^^{x;z)-[fi^u]. (3) 



The self duality of F^^ implies D"^ D = —aoD^D^. Accordingly, {DlDz)~^{x,x') commutes 
with all cr^, from which the self duality of F follows. The 'original' gauge potential, A^{x), 
can be recovered by Nahm transforming A^{z) 

A^^ix) = I d'z r\x-, ^)^^r{z; x), (4) 

where the ip^lz-jx), p = 1,2,. ..,N are an orthonormal set of zero modes of -Dj.(A) = 
—crjj^Dt^^A) with Dt^ = d/dZfj, + — ix'^. In other words the four torus Nahm trans- 

form is involutive. 

Formally, one can obtain the x Nahm transform by taking two of the periods, say 
Lq and L3, to be infinite. Given an SU (N) instanton periodic with respect to Xi —>■ Xi + Li, 
X2 ^ X2 + L2 its Nahm transform is 



= / d'x^|J'\x■,z)^^P^ix■,z), /x = l,2, 

Af{z) = / d^xilj'\x;z)iXf,^\x;z), /i = 0,3, (5) 

where the z) with i = 1, k are orthonormal zero modes of Dl{A). Note that we may 
gauge Zq and z^ to zero. By assumption the zero modes z) are square integrable. This 
does not guarantee that the integrals in (0) exist. However, the integrals only diverge at 
special z^ values where the zero modes do not decay exponentially. The number of these 
points turns out to be N, i.e. it is determined by the gauge group of the instanton. To 

summarise, the Nahm transform maps doubly periodic SU (N) instantons with topological 

~ 2 

charge k to self dual U{k) potentials on the dual torus, T , with singularities. 

The dimensionally reduced self-duality equations (or Hitchin equations) take a particu- 
larly simple form in complex coordinates 

y = zi+iz2, y = zi-iz2, (6) 

with derivatives dy = \{dzi — idzi), dy = + idz2)- Combine Ai{z) and A2{z) into a 

complex gauge potential. Ay = ^{Ai{z) — iA2{z)) and Ay = ^{Ai{z) +iA2{z)), and form a 
'Higgs' field out of the remaining components, = ^{Ao{z) — iAs{z)). The self-duality 
equations read 

[Dy,^=0, Fyy=[Dy,Dy] = [<^,<!>^]. (7) 

The reader should be aware that we will occasionally use both the cartesian z and complex y 
coordinates in the same equation. The Nahm transform (jSJ is a straightforward limit of the 
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version. The inverse transform is a different matter. That is, given a (singular) solution of 

~ 2 

the Hitchin equations on the dual torus, T , how does one recover the corresponding doubly 
periodic instanton? To our knowledge there is not a mathematically rigorous treatment of 
this issue in the literature. Our approach will be to study specific simple solutions of the 
Hitchin equations and investigate whether the number of normalisable zero modes of D^A) 
matches the number of singularities of A{z). When this is the case we can then ask whether 
the natural analogue of (@]), namely 

^r(^) = j_/zr\z-x)^^r{z;x), (8) 
provides doubly periodic instantons. 

3 Charge One 

In the one-instanton case the Nahm potential is abelian and the Hitchin equations reduce to 

Fyy = 0, dy<l> = 0, (9) 

SO that $ is anti-holomorphic in y. As we shall see. Ay can be taken as holomorphic. Since Ay 
is periodic it is an elliptic function. It is well known that they have at least two singularities. 
We expect the number of singularities to correspond to A^. Let us stick to the simplest case 
N = 2. Thus Ay will have two simple poles ^. Consider the ansatz 

Ay = dy(f), Ay = -OyCf), {W) 

which gives Fyy = —2dydy(l). It follows from (jH)) that (f) must be harmonic except at two 
singularities. A suitable satisfies 

{di + dl) Hz) = -2nK [6'{z -uj)-6'{z + u;)] , (11) 

where k is a constant and ±uj are the positions of the two singularities (we have used 
translational invariance to shift the 'centre of gravity' of the singularities to the origin). 
The delta functions should be read as periodic (with respect to zi zi + 2tt/Li and 
Z2 ^ Z2 + 2t[/L2). Physically, the Nahm potential describes two Aharonov-Bohm fluxes of 
strength n and — k threading the dual torus. They must have equal and opposite strength 
to ensure a periodic A. We may assume that k lies between and 1 since it is possible via a 
"'^The Weyl equations for Nahm potentials with higher poles do not admit normalisable zero modes. 
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(singular) gauge transformation to shift n by an integer amount (under such a transformation 

~ 2 

the total flux through T remains zero). One can write explicitly in terms of Jacobi theta 
functions (see for example |rij) 



K 

2 



log 



^((?/+^i+^^2)i^ + i + it,t) 



v 



\ 

2 H z — yy -y) - ^^2Li 



TT 



, (12) 



where y and y are the complex coordinates introduced in 
as 



The theta function is defined 



9{w, 



T 



iirri'^ T+2iTinw 



Imr > 0, 



(13) 



and has the periodicity properties 6{w + 1, r) = 9{w,t) and e{w + r, r) = e-'''^-^^''"0{w, t) 
In each cell 9{w,t) has a single zero located at the centre of the torus (ly = | + |r). We 
have chosen the constant term in (fT^ so that the integral of (p over the dual torus is zero. 
This renders (j){z) an odd function, (p{—z) = —4>{z). It is easy to see that inserting (|T^ into 
()1U|) yields a holomorphic Ay with simple poles at the fluxes. 

What about the Higgs field? Since we are aiming for an SU{2) instanton it must, like 
A{z), have two singularities. The singularities must be at the same two positions as in Ay or 
else we have a total of four singularities which means we are considering a (rather special) 
SU{4) instanton. The simplest way to arrange for Ay and the Higgs to have singularities at 
the same points is to choose them to be proportional 



= ad, 



(14) 



where a is a complex constant. Not all possibilities are exhausted, since while the poles 
must coincide the zeroes need not. This remaining ambiguity corresponds to the freedom to 
add to $ a complex constant. However when we insert our Nahm potential into the Weyl 
equation such a shift is equivalent to a translation of xq and X3. This Nahm potential was 
derived in [Hj via the ADHM construction. Here the parameters k, and a are related to the 
'size', A of an instanton centred at x. 



0, 



ttA^ 
L1L2 



(15) 



Although we do not directly use the ADHM construction in the present paper this relation 
proves useful in interpreting our results. The Nahm potential we have described has five 
parameters; Ui, UJ2 fixing the flux separations, n the flux strength and the complex parameter 
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a specifying the Higgs field. Then, with the four translations in x M^, the corresponding 
doubly periodic instanton has a total of nine parameters. 

To reconstruct the charge one S\J (2) instantons corresponding to our simple Nahm po- 
tential we require the two zero modes of the Weyl operator 

- -Dl{A) = ( ^^^^ i(^^~^dy^ dy + dy(P - 
^ I 9y — dy(j) — |x_L — iaK~^dy(f) 

~ 2 

In addition to the complex coordinates y, y on T we have introduced two sets of complex 
coordinates for T X M^; in the 'parallel' directions x\i = Xi + ix2^ xn = Xi — ix2, and in the 
'transverse' non-compact directions x± = Xq + ix^, x± = Xq — ix^. In this paper we shall 
concentrate on the special case a = 0, i.e. zero Higgs. This means that the corresponding 
instanton will be radially symmetric: the action density depends on xi, Xi and r = only. 
These radially symmetric solutions are a seven parameter subset of the doubly periodic SU (2) 
one-instantons. 




4 Nahm Zero Modes 

From now on we stick to the zero Higgs case (i.e. a = 0). Even with this specialisation 
the Weyl equations are still rather forbidding. However, there is one special case which is 
immediately tractable: when xj, = the Weyl equation decouples and the two zero modes 
have a simple form jHj 

i,\z-x)=\ _ ° I, ^\z-x)= \ ^ ' 1 (17) 




where G±{^z^ are periodic Green's functions satisfying 

{-%dy - G^{z) = l6\z), {-tdy - G4z) = l6\z). (18) 

Here we have used both the z and y coordinates in the same equation; this possibly confusing 
usage will continue until the end of section |S| G- {z) has a theta function representation 

"a i+ife.t)' 

with 9'{w^ t) = dyj6{w, r). The corresponding result for G^{z) can be obtained via G^{z) = 
G*_{—z). The zero modes have square-integrable singularities at both fluxes: near z ^ u 
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we have l^J^]"^ oc \y ~ cui — iuj2\'^^'^~^^ and l^/'^P oc \y — cui — iuj2\~'^'^. Near the other flux the 
singularity profiles of the two modes are exchanged in that oc \y + Ui + iuj2\~'^'^ and 

oc \y + uoi + iuj2\^^'^~^^ ■ By virtue of the second theta function in the numerator of (|T9|) 

~ 2 ...... 

il)^{z] x) has a single zero in T at y = —ui — iuj2 — 2Tiix\\/ {L1L2) (see also T^). Similarly, 
'ip^{z] x) has a zero aX y = Ui + iuj2 — 27rix||/ (L1L2). 

5 Soft Zero Modes 

When x\\ = the Green's functions G± do not exist. This does not mean there are no zero 
modes, indeed one can see that 

^l;\z;x = 0)= I ^p^{z;x = 0)=\ ^ 1, (20) 

are solutions of the X|| = Weyl equation. has the expected square- integrable singularity 
at 2; = —uj, but for z = uj, iIj^{z]x = 0) is zero. On the other hand tlj'^{z,x = 0) diverges 
at z = uj but not at z = —uo. By choosing x\\ = we have effectively moved the zeroes 
of the zero modes to flux locations. There is another way of doing this; choosing x\\ = 
—iLiL2{uJi + 1002) /it brings the zero of ip'^ to the other flux (and similarly for ip"^). Thus 
there are two x\\ values where each Nahm zero mode diverges at only one flux. We wish to 
investigate whether these 'soft' points can be interpreted as the locations of some kind of 
constituent. 

The singularity profiles of ip^ and ip"^ are exchanged under the replacement k 1 — k 
suggesting that the constituents are exchanged under this mapping. That is, if there are 
indeed lumps at the two points, then n ^ 1 — n swaps the two lumps. The following result 
formalises this idea 

F^^{x\\,x^, k) = V~'^{x)F^^ ^-a;|| + ^^T^(t^i + ^^^2), -x±, 1 - V{x), (21) 

where V{x) is some U{2) gauge transformation. The proof of Eq. (PT|) goes as follows: make 
the change of variables z — > —z in (jHJ. The zero mode ipP{—z;x) satisfies the same Weyl 
equation as ip^^z; x) except that the signs of k and the x^ are flipped. However, —n does 
not lie between and 1. Under periodic gauge transformations A{z, —k) and A{z, 1 — k) are 
equivalent up to a constant potential 

A^iz, -k) = U-\z) (d^ + A^{z, 1 - k) + 5^) U{z), (22) 



where ttBi = ~iLiL2UJ2 and 7ri?2 = iLiL2UJi. Explicitly 

U{z) = e^-'^L,L2{y-y)/n ^23) 



X 



In the Weyl equation can be absorbed into xi and X2- Thus U ^{z)iIjP{z, —x\\ —iLiL2{uJi + 
iuj2)/T^, — satisfies the same Weyl equation as ip^{—z;x\\,x^). 

6 Green's Function Approach 

For the transform, the original field strength can be written 

Fll{x) = [ d'z [ d'z'V\z;x)a, iDlD,)-\z,z') alr{z-,x) - [/i ^ v], (24) 

and as with the corresponding expression for F^{,{z), equation Q, the self-duality of FJ^^^x) 
is manifest because the Green's function {DlD.j.)^^{z, z') commutes with all cr^. Starting 
directly from the inverse Nahm transform (0)) we have 

K'i^) = i d^zd,r\z-x)d,ilj\z-x) (25) 
+ / d'^zi)P\z]x)d^i)''{z]x) I d^z'ij'^z';x)d^ij''{z';x)-[fi^u], 

where 9^ = d/dx^. To get from (j25j) to (j24p we require the derivative of the Nahm zero 
modes with respect to the coordinates of 

A_^i = -^D,{A) (^Dl{A)D,{A)y\lr + r R^i^). (26) 

for some N x N matrix valued vector field on the four torus, R1[{x). This clearly satisfies 
the equation obtained by differentiating the Weyl equation, DI{A)iIj'^{z; x) = 0, with respect 
to x^. Multiplying from the left with tp^"^ and integrating over the dual torus yields 

RP'^ix) = AP\x), (27) 

which together with equations and give (j^ . 

Our analysis of doubly periodic instantons has been based on the Weyl equation and 
the inverse Nahm transform (jS)). These are the exact counterparts of four torus equations. 
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Therefore one might expect all the results quoted above to carry over directly to x M^. 
There are, however, a number of important differences between the two cases. With regard 
to the X construction, the Green's function {Dl.{A)D,j.{A))~^{z, z') does not commute 
with the 0"^ matrices. Away from the singularities of A{z), Dl.{A)D^{A) commutes with 
the (T^. At the singularities it has source terms which are not proportional to ctq. For the 
f/(l) Nahm potential considered in section 3 (as usual we set a = 0) 

DlD, = -ao{D^ f + iTiiKo^ {b''{z -uj)-b''{z\ cu)] , (28) 

indicating that {D^D^)^^ does not commute with a\ and o^- 

We will show that the naive X limit of ^ does not apply, so that {D\D:,)-^ 
not being proportional to ctq does not signal a breakdown of the inverse Nahm transform. 
Starting with the inverse transform (jH}, the field strength can be written in exactly the same 

~ 2 ~ 4 

way as in except that the integral is over T rather than T . The derivative formula 
(|2o|) also applies without modification to the T zero modes. What does not carry through 
is equation (j27|) . Much as in the four torus analysis i?^ can be determined by multiplying 
the X analogue of ()26|) by ^^^^ and integrating over T . 

Al\x) = -ij^r^DM) [Dl{A)DM)Y'<ylr + Rl'. (29) 

Naively one would expect the integral on the right hand side to vanish as ip^'^ is a (left) 
zero mode of D^{A). However, this property only leads to vanishing i^'^'^ Dx{A)(j) for suffi- 
ciently smooth functions (j). In (|29|) ip^'^ and the Green's function {Dl.Dr^)~^ have coincident 
singularities. 

Retaining the integral yields a modified field strength formula 

F',l{^) = i d'z f dh'rHz;x)a,fiz,z';x)alriz';x) - [f, ^ z/], (30) 

with 

f{z,z';x) = {DiD,y\z,z') (31) 
- /_ d's{DlD,y\z,s)Dl{A)r{s;x) 

X / dh'rHs';x)D,{A) {DlD,)-\s',z'). 

This object does commute with the o"^ consistent with self-duality. The field strength for- 
mula ()30|) is similar to those one can derive for instantons and calorons via the ADHM 
formalism. 
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To make these considerations more concrete we see how they apply to our simple U{1] 
Nahm potential. In this case 



(32) 



^ * l\X±\^ + {-dy + dy(f) + ^X\\){dy + dy(t) - ^X\\) 

When X this is invertible. The inverse has the form 

2 {DlD^y' iz,z') = ((To+«(T3)e-*(^)ir+(^,z';x)e-*(^') + (ao-2a3)e*(^)ir„(^,^';x)e<^(^'), (33) 

where the K±{z, z'; x) are finite for all z and z' (unless x\\ = x± = 0). Because of the 
exponentials {DlDx)~^{z, z') is singular if z or z' equals ±u!. In the x± = slice the K± 
have simple integral representations 

K±{z,z';x^ = 0) = [ dhG±{z - s)e^^'t'^'^G^{s - z'), (34) 

where G± are the Green's functions defined in section 4. The x± = zero modes of section 
4 can be rewritten as follows 

/ e'^^'^K(z,u;x) , 
ij\z;x) = -D^iA) I ' ^^^^ 



e 



^'^K+{z,-oj]x) 



In fact, these are also valid for x±_ ^ 0. The modes are orthogonal and can be normalised 
by dividing them by ^/p, where 

p{x) = Kj^{—ui,—uj]x) = K_{(jj,uj]x). (36) 

Inserting the zero modes into ()31|) enables us to prove that / commutes with 0"i,cr2- This 
is detailed in appendix 1X1 One can also show that f{z, z'\ x), unlike {D\D^~^{^z^ z'), is well 
behaved at x = 0. As 2; or 2;' approaches ±uj, f tends to zero. 

We have seen that it is possible to construct the Nahm zero modes out of objects, K±, 
entering the inverse of DlD^. The gauge potential can be computed from these zero modes. 
The final result is very simple 

T3 ^* 
= -—9^11 log p - 27ri(ri - ir2)/tp9j;^ — , (37) 
Z p 

T3 ^* 
= -—d^^\ogp + 2ni{Ti-iT2)Kpds^^ — , 
Z p 

11 



where 

^{x) = K_{lj, —uj; x), i'*{x) = K_{—uj,Lj;x). (38) 

Note that p is dimensionless, real and periodic (with respect to Xi — » Xi + Li, X2 — > X2 + L2), 

while u is dimensionless, complex and periodic up to constant phases. Eq. (j37j) can be 

checked component- wise (see appendix IH]) . 

We have given closed forms for the K± in the special case x± = 0. An explicit construct 

valid beyond this two-dimensional slice would immediately provide the exact Nahm zero 

modes and gauge potential. Note that the Green's functions K± are radially symmetric. 

Therefore the functions p and u are as well. They can be expressed as power series 

00 00 
P = $^Pnkx|'", = J]z/„|xxP". (39) 

n=0 n=0 

7 Field Strength Components 

In section 5 we saw that the Nahm transform singled out two points in the x± = slice. This 
suggests that they may be the locations of some kind of constituent. As they are isolated 
points in x rather than lines it appears that we are not dealing with monopole con- 
stituents. Indeed if we take (fT3j) at face value the constituents appear to be BPST instantons, 



one with size A = ^/kLiLqJt: at X|| = x_l = 0, the other with size A = a/(1 — k)LiL2/'k at 
the second soft point. This description shares some features with early attempts to describe 
SU{2) instantons in terms of 2k 'instanton quarks' [TB] . 

To see if this is realised we must compute the field strengths. This can be done explicitly 
in the a;_L = slice. Inserting the power series (j39|) into the ansatz (|37|) it is straightforward 
to compute F^^^s^^ and 



F - 

^ii^ii 



Tsd^ .dsAogpo, F^ s^ = 2'Ki{ri + iT2) Kpod^ —. (40) 



=0 ' ■ ■ "po 



The other components are fixed by self-duality, i.e. + ^x^x^ = and F^^^^^ = 0. Both 

Po and z/q diverge at x\\ = 0, but the field strengths should not. A careful analysis of po 
and uq in the neighbourhood of x\\ = shows that the field strengths are well defined at 
this point. Alternatively, one can just note that p and u are well behaved at the second soft 
point and invoke (j^D). Using these results the action density can be computed 

2 



iTr F^,F'^n.^=o = 16 (s,,,^^,, logpo)' + 1287r\Vo 



52 ^ 

^"Po 



(41) 
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From ()34|1 one can derive integral representations of po and uq where the integrands are 
expressible in terms of standard functions. These representations allow us to plot the action 
density within the x_l = slice for different values of the parameters k, ui and 002- The two 
periods Li and L2 can also be varied. However, the physically important quantity is the 
ratio of the lengths 

a = ^, (42) 
since scaling identically both lengths equates to a trivial scaling of the action density: 

p^TiF^Px; pLi,PL2; p-^u^, p-^cu^; k)=TiF^{x; L^, uo^.lo^, k). (43) 

In jHj plots for different n in the equal length case, a = 1, were presented. There we found 
that the action density has either one or two peaks. A single peak is observed when the two 
soft points are quite close or when one size is somewhat larger than the other. In the latter 
case the smaller sized instanton dominates. Even when two peaks are resolved there is quite 
a strong overlap. In fact such an overlap is 'necessary' since otherwise the instanton would 
have charge two rather than one. The constituent locations then correspond to the 'cores' 
of the overlapping instantons. 

Here we also consider rectangular tori with a 7^ 1. Without loss of generality we set 
1/2 = 1 and Li = a. To begin with consider k = | where the two constituents are identical. 
Now choose (001,002) = {^na^^ , ^n) . This means that there is a core at the centre of the 
torus and another of equal size at the corners. Starting with a = 1 increase the value of a 
thereby stretching the torus in the xi direction. 

In fig. ^ the action density is plotted for some values of a ranging from 1 to 3. The first 
plot (a = 1) is just the equal length case clearly showing two instanton- like peaks exactly at 
the expected locations. The second plot (a = |) shows that as the torus is stretched the two 
cores become stretched as well; they no longer have the approximate spherical symmetry 
of the a = 1 case. On stretching the torus further the now deformed cores seek to overlap 
with the periodic copies of themselves. This has already started to happen in the third plot 
where a = 2; the soft point locations are no longer maxima of the action density. Finally, 
at a = 3 the self-overlap has generated two monopole like objects; the solution has a very 
weak dependence on X2 and there are peaks at two values of Xi. Note that the two monopole 
'worldlines' do not cross the soft points. Increasing a more weakens further the dependence 
on X2- In the limit a ^ 00 a periodic monopole (i.e. a self dual solution on x M^) will 
form. However, this is a singular limit. 
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Figure 1: action density — ^Tr for x± = 0, k = ^ and ui = \tx, uj2 = ^na for a = 1, |, 2, 3. 

What is happening is that the doubly periodic instantons are becoming more caloron- 
hke (caloron-hke objects have been observed recently in lattice-based approaches ITH] 
). It is helpful to recall some of the basic features of the simplest caloron solutions. An 
SU(2) caloron with charge one comprises two monopole constituents separated by a distance 
7rA^/L2 (taking X2 to be the compact coordinate). Here UJ2 (there is no ui) determines the 
mass ratio of the two monopoles. For the x problem ttA^ = HL1L2 suggesting that the 
separation of the monopoles should be simply uLi for k < ^ and (1 — k) Li for k > i. In 
fig- Q we took n = \ and the worldlines are indeed separated by a half-period. In fig. |21 we 
have taken u as above but increased to |. The first plot is the equal length case showing 
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a peak at the centre of the torus. At the origin of the torus, the location of the larger 
constituent, there is a broader but much shallower peak. The second plot is for a = 3 which 
clearly shows two monopoles. Moreover, the distance between them is consistent with 
Though the dependence on X2 is not quite as weak as for the corresponding plot of fig. [TJ 
away from k = | more stretching is required to access the monopole regime. 




Figure 2: action density — ^Tr for x± = 0, n = ^ and ui = ^n, uj2 = ^vra with a = 1,3. 

We have argued that for a at or near to 1 the constituent core picture best describes 
the action density while for larger (or smaller) a almost static monopole-like objects form. 
In that respect, the instanton core and monopole constituent pictures are complementary. 
However, we have an example displaying features of both descriptions: 




Figure 3: Plot of action density for = 0, k = |, a = 1 and uii = ^n, U2 = 0. 
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Above we plot the action density for k = ^, Ui = and 002 = 0. This configuration 
shows essentially no dependence on X2 even though the lengths are equal. What seems 
to have occurred is that the two cores at x\\ = and x\\ = |i have merged into a single 
monopole. Unlike in the previous examples the monopole worldline passes through the soft 
points. That there is only one worldline is not inconsistent with the constituent monopole 
picture since by virtue of ti;2 = the 'second' monopole should be massless. The monopoles 
observed in figs. ^ and |21 have equal masses in accord with the symmetric value of 002- 

8 Charge One Half 

When K, = \ the two instanton cores are identical. This has an interesting consequence. If 
we choose the constituent locations so that they are separated by half periods the charge 
one instanton can be 'cut' to yield two copies of a charge | instanton (see also [H]). This 
happens when {lji,lj2) is (|7r/Li,0), (0, ^tt/L2) or {^n/Li, |7r/L2), see fig. IH After cutting 
we have a twist Z12 = — 11 in the half torus. But to produce a half integer topological charge 
another (non-orthogonal) twist is required; this is most simply achieved with Zq^ = —11. 
Such a twist would have the novel feature of being associated with the non-compact Xq and 
X3 directions. Far away from x± = the potential must be pure gauge 

A^ix) ~ V-\xu X2, e)d^Vixi, X2, 9), r ^ 00, (44) 

where x± = re*^. The non-compact twist translates into a double valued gauge function, 

V{Xi, X2, e + 2TT) = -V{Xi, X2, 9). 

The examples considered in fig. Q are all 'doubled' half instantons. Here the cut tori are 
diamond shaped (see fig. 4) apart from the first (a = 1) plot, where it is square. This case 
has another interesting feature. Compare the constituent locations with the half instanton 
obtained via Li = \/2, L2 = l/V^, lui = 0, uJ2 = vr/y^. On cutting both yield charge | 
instantons on square tori with length By analogy with the four torus case, for a given 

set of twists and periods, charge | instantons are expected to be unique up to translations. 
Yet these two cases have different values of a, one has a = 1 where the core picture should 
work while the other (a = 2) is expected to have some monopole characteristics. As can 
clearly be seen from the first plot in fig. ^ the action density has maxima at the soft points. 
We have computed the a = 2 case and found maxima not coincident with the soft points. The 
soft points correspond to saddle points of the action density. But allowing for translations 
the two half instantons are identical. What seems to be happening is that when a = 2 the 
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Figure 4: When k — the choices uji = ^tt/Li, uj2 = and ui = ^n/Li, U2 = \t^IL2 
yield constituent separations allowing the torus to be cut (dashed line) to yield charge | 
instantons. 



monopoles are still not developed, and in this intermediate state there are lumps of action 
density which while not at the soft points are identical to those at the soft points of the 
a = 1 case. More precisely the two action densities satisfy 

l^Fl^^{x,,X2) = IVF„li(^(xi + - i - X,) - \). (45) 

The a — 1 plot is obtained by a 45 degree rotation of the a — 2 plot followed by a translation. 
This reflects the fact that the two Nahm potentials are also related by a 45 degree rotation. 

9 Asymptotics 

In this section we consider the large \x^\ limit of the zero mode equations, and deduce the 
large \x±\ behaviour of the gauge potential and field strength. This allows us to check that 
the charge | instantons have the expected properties. In particular, we see that the gauge 
potential has asymptotic behaviour commensurate with a twist, Zq^ = — H. Furthermore, 
the exponential decay of the action density is derived. In fact, all the — \ solutions have 
both the twist and exponential decay even if they are not 'doubled' | instantons. 
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When |a;_L| is much larger than the two periods the zero modes of D^A) become strongly 
localised about the two flux singularities. In this regime we can approximate the zero modes 
with solutions. A single flux in corresponds to = —^nlogyy. The Weyl operator is 
then 

When x± 0, Dl has a normalisable zero mode. In the special cases k = ±| the normalised 
modes take rather simple forms: 

^pi(z) = — , 1 , i(z) = — , ^ , (47) 
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X±\ 



y 

where x ■ z = xiZi + X2Z2. For sufficiently large |a;_L| we can take 

iP^z) =^{z- uj), ij\z) = ij_.{z + u), (48) 

and extend the integration region from T to M^. For this approximation to work it is 
important that the flux separation 2\uj\ is less then both dual periods, 2tt/Li and 27r/L2. 
This is because we assume that the 'interference' between the fluxes dominates the effect of 
the periodicity for large |a;_L|. Note that 

d^zilj\{z-uj)ilj_i{z + uj) = 0. (49) 

ij2 2 2 

Inserting the zero modes (j^Hj) into the inverse Nahm transform (jHJ produces a U{2) potential; 
we discard the pure gauge U{1) part. The remaining SU{2) piece turns out to be in a different 
gauge to that implicit in section IHl An explicit computation (see appendix C) shows that 
the SU (2) potential has asymptotics corresponding to a twist and the action density has the 
exponential decay 

- \W,^F,^ ^ ^,-4|.||.xl (50) 

We can apply this to the charge | case. Consider such an instanton in a square box of length 
L. This is a one instanton in the doubled box with Li = L, L2 = 2L. Here ui = 7r/(2L), 
UJ2 = 0, and so the | instanton has the fall-off 

- ^TrF,.F,. ~ ^e--l-l/-. (51) 
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For a rectangular box simply replace L with the longest length. In j2! a decay exponent of 
6.5 was reported for a unit box. This is within 5 percent of the analytic value 27r. The decay 
formula is reminiscent of similar results for Bogomolnyi vortices in the abelian Higgs model. 
There the decay exponent is straightforward to obtain while an analytic computation of the 
prefactor is still an open problem (see however ITH]). 

One can repeat the analysis for k 7^ |. The corresponding zero modes are more 
complicated; they can be written in terms of modified Bessel functions. Here the asymptotics 
of the gauge potential does not correspond to a four torus twist. More precisely, p has the 
power decay 

p = — + exponentially decaying remainder, (52) 

l^^-L I 

where C is a (dimensionful) constant. This translates into a twist only if k = |. 

The question of the form of the decay was an important consideration in T]. There 
it was assumed that doubly periodic instantons have quadratically decaying field strength 
components. This translates into a quartic decay of the action density. Yet the solutions 
we have discussed show a much more rapid exponential decay. This is not necessarily a 
contradiction since our discussion of the asymptotics was confined to the radial solutions 
whereas was concerned with the whole charge one sector. The more generic non-radial 
solutions may decay quartically. Indeed the two monopole constituents are expected to be 
separated in the non-compact directions. This can give rise to an algebraic tail in the action 
density. 

10 Outlook 

We have developed an instanton core picture of the simplest doubly periodic instantons. 
Numerical calculations of the action density correlate well with this picture for square tori. 
However, in rectangular tori the solutions become caloron-like when the ratio of the two pe- 
riods is increased; the cores merge with periodic copies of themselves in the short direction to 
form monopole-like tubes of action density. The basic properties of these monopoles, such as 
their mass ratios and spatial separation, follow the pattern of the charge one SU{2) calorons. 
It would be interesting to develop this monopole constituent description in a more direct 
manner. One approach would be to exploit the fact that the radially symmetric solutions 
considered here fall into a class of axially-symmetric multi-calorons discussed recently JHl if 
one allows for infinite topological charge. 
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There are a number of obvious ways to extend the results in this paper. An exphcit 
treatment of the 7^ zero mode equation is still lacking. In the non-radially symmetric 
case (a 7^ 0) the zero mode equations have not yielded solutions for any points in x M^. 
Even in the absence of explicit zero modes it might be possible to extract some information 
about possible instanton core or monopole constituents. A more straightforward extension 
would be to generalise the results regarding radially symmetric one instantons to SU{N). 
That these solutions decay exponentially indicates they can be compactified further to T"' x M 
or even T^. The latter is only possible for the where the asymptotics of the gauge 

field correspond to a twist. 
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A Properties of / 

Inserting the zero modes ()35p into the definition of / (j^T|) yields 

f{z, z'- x) = i(^^o + ia,)e-^^-^^g+{z, z'- x)e-<^(^') + i(^o - ia^)e'i'^-^^ g.{z, z'- x)e'^^''\ (53) 
where 

9±{.z, z'] x) = K±{z, z'; x) —K±{z, =Fu;; x)K±{tuj, z'; x). (54) 

p{x) 

It is not obvious that / commutes with ai and a2- This only holds if the coefficients of the 
projectors ± ia^ match, that is if 

e-^^'^g+iz, z'- x)e"'^("') = e'^^''>g^{z, z'- x)e'^'^''\ (55) 

A proof of ()55|1 for x_l = was given in 8J . It is sufficient to show that the left and right hand 
sides of ()55|) have the same asymptotics at the fluxes since when z, z' ^ -^uj they satisfy the 
same differential equation. Consider g±{z,z';x) in the neighbourhood ol z = uj. It follows 
immediately from (j54|) that g_{ijj,z'\x) = whereas g-^^{uj, z'; x) is non-zero. This does not 
contradict ()55|) since e'^^^^ diverges (oc \y — cui — iuj2\^'^) at z = cu. A short calculation gives 

g. (., z;x)r. K_ (u;, z;x). (56) 
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This is compatible with if 

g+{uj, z'; x) = —K^{uj, z'; x). (57) 

A'n-Kp[x) 

To check this, note that for z' ^ ±uj both sides satisfy the same differential equation and then 
examine them in the neighbourhoods of z' = ±uj. From this discussion of the asymptotics of 
g± it should be clear that f{z, z'\ x), unlike {D\Dx)~^{z^ z') is always finite. Indeed /(z, z'\ x) 
tends to zero as z or z' approaches a flux point. 

B Computation of ( IHTj ) 

Here we outline the computation of Ax^^\ the other components can be dealt with in much 
the same way. Inserting the zero modes into the inverse Nahm transform yields 

= ^l_^d'zK.{u;,z-x)e'^''^'> (^DUA)d.^^D,iA)y' ^K.{z,u- x). (58) 
Using (L)ta^||D)ii = d^^^iDWy^ and the representation (33) of (DW)-^ 

K = ^ / d'zK_{u,z;x)d,^^6\z -u)^ = -^d.^^logp, (59) 
as required. For A^^^ we proceed similarly 

All = ^ dhK+i-u;,z;x)e-'^^^UDliA)d,MA)) ^K4z,uj;x) (60) 



= i — / d z K+{—u!^ z; x)K_{z,uj; x). 
P Jf^ 

The key step is to use ^7\\ 

Al]=47riKX^ dhK+{-uj,z]x)e-^'''^'^g+{z,uj;x). (61) 

I Jf2 

Now use (j54j) to write Qj^ in terms of and the formula 

/ d^sK+iz, s; x)e-2*Wir+(s, z'; x) = --^K+{z, z'; x), (62) 
to get the result Af- = — 47rm9j.^(z/*/p). 
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C Exponential Decay of the Action Density 

We consider the A^j.^ derived from the asymptotic zero modes given in section 9. A straight- 
forward calculation gives 

K - < = 0. (63) 
The off-diagonal components are more involved: 

Al = I ^.iU^ + - = ^e-"(/ -H J). (64) 



Here 



^2^(M^(^^^^)-ie-l"^l(v^+^^) (65) 



i?2 2/2 



and 



J=- [ rf2^e-l^^l(v^+v^) (66) 



X 



. _ s_i (2/12/1)4 ,(2/22/2)4, _ . 1 I (2/12/1; 
(2/22/2) 4 — z \x±\ — z (2/12/1)* - F± I (2/22/2) 



2/1 2/2 2/1 



where yi = y — ui — iuj2 and 1/2 = y + ^1 + i^2- As is increased the integral / is 
dominated by a small neighbourhood about the line-segment joining the two fluxes. For 
example, if u;2 = we have e'l^^K^^^+v^^) ^ 7r2 |x_Lr2 |u;i|"2(u;2 _ ^2^2e-2|^±ll^il5(2-^) for 

—uoi < zi < uoi. Integrating along the line-segment gives 

j^^![WWie-H|.J. (67) 
UJl — iuj2 

J exhibits a similar exponential decay. is simply x ±e^'^^'^ J / {2'k\x ±\) , and so 



1 

1 2ix-up._ (^T^ ,2»x...,-2M|x|| (^l +^^2)|t^| 2 



^^l^^l 7r2|x±|2 

We also have 

F^l., = O(e-^l-ll^'-l). (69) 
Now compare with the corresponding field strengths generated by the ansatz ()37p in sectionlHl 

F^l,^ = 2mpd,^^d,^j, Fl\,^ = d^^d,^ logp + (27rp)29.,|^5,,|^. (70) 



Consider 

C . 1 



P = 1 77. ^* = -vr-2e"*"-"e-'l"ll"^Vr^|2;±r^ (71) 
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where t is some real number and C is a constant. Neglecting sub-leading terms the ex- 
pressions for Fx^x^ agree up to a phase, \/x±/x± = e*^, which can be compensated for by 
the double- valued gauge transformation V = e5*^'^3_ Under this gauge transformation the 
diagonal part of A^^ becomes non-trivial 

= -^ + 0(e-^l"ll"^l). (72) 

_|_ 

which agrees with (jTTj) if t = 1 in (ffT|) . This is exactly what is required to generate the twist 
= —11. Using the large form of p and u it is straightforward to obtain the decay 
formula (fHOjl . 
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